arXiv: 1502.04711 v2 [astro-ph.GA] 18 Feb 2015 


Astronomy Astrophysics manuscript no. OrionGEMS ©ESO 2015 

February 19, 2015 


The Orion Fingers: Near-IR Adaptive Optics Imaging of an 

Explosive Protostellar Outflow 

John Bally 1 , Adam Ginsburg 2 , Devin Silvia 3 , and Allison Youngblood 1 

1 Department of Astrophysical and Planetary Sciences, 

University of Colorado, UCB 389, 

Boulder, CO 80309 

e-mail: John.BallyQcolorado.edu 

2 ESO Headquarters 

Karl- Schwarzs child- Str. 2 

85748 Garching bei Munchen, Germany 

3 NSF Astronomy and Astrophysics Postdoctoral Fellow, 

Department of Physics and Astronomy, 

Michigan State University, 

East Lansing, MI 48824 


February 19, 2015 


ABSTRACT 

Aims. Adaptive optics images are used to test the hypothesis that the explosive BN/KL outflow from the Orion OMC1 
cloud core was powered by the dynamical decay of a non-hierarchical system of massive stars. 

Methods. Narrow-band H 2 , [Fe n], and broad-band K s obtained with the Gemini South multi-conjugate adaptive optics 
(AO) system GeMS and near-infrared imager GSAOI are presented. The images reach resolutions of 0.08 to 0.10 // , close 
to the 0.07" diffraction limit of the 8-meter telescope at 2.12 pm. Comparison with previous AO-assisted observations 
of sub-fields and other ground-based observations enable measurements of proper motions and the investigation of 
morphological changes in H 2 and [Fe n] features with unprecedented precision. The images are compared with numerical 
simulations of compact, high-density clumps moving ~ 10 3 times their own diameter through a lower density medium 
at Mach 10 3 . 

Results. Several sub-arcsecond H 2 features and many [Fe n] ‘fingertips’ on the projected outskirts of the flow show 
proper motions of ~300 km s -1 . High-velocity, sub-arcsecond H 2 knots (‘bullets’) are seen as far as 140" from their 
suspected ejection site. If these knots propagated through the dense Orion A cloud, their survival sets a lower bound 
on their densities of order 10 7 cm -3 , consistent with an origin within a few au of a massive star and accelerated by a 
final multi-body dynamic encounter that ejected the BN object and radio source I from OMC1 about 500 years ago. 
Conclusions. Over 120 high-velocity bow-shocks propagating in nearly all directions from the OMC1 cloud core provide 
evidence for an explosive origin for the BN/KL outflow triggered by the dynamic decay of a non-hierarchical system of 
massive stars. Such events may be linked to the origin of runaway, massive stars. 
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1. Introduction 

The BN/KL region behind th e Orion Nebula, loca ted at 


a distance of about 414 pc (Menten et al. 12007), con¬ 
tains a spectacular, wide opening-angle, arcminute-scale 
outflow emerging from the OMC1 cloud core. The flow 
is traced by the millimeter and sub-millimeter emission 
lines of molecules such as CO, CS, SO, SO 2 , and HCN 
that exhibit broad (> 100 km s" 1 ) emission line wings 


(Kwan & Scoville 1976 

Wiseman & Ho 1996 Furuya & 

Shinnaga 

2009), high-vel 

ocib 

y OH, H 2 O, and SiO maser 

emission ( 

Cenzel et al. 1981 

| Greenhill et al. 

TT J LS 2 . 

1998), and 

l* 1 0 a 


pm 

Fe ill emission ( 

'Allen & Burton 1993[ Colgan et al. 

2007 

Nissen et al. 2( 

307| |Lee & Burton 2000; Bally et al. 

2011 

). The OMC1 outflow has a southeast (red-shifted) to 


northwest (blue-shifted) axis and contains at least 8 Mq of 
accelerated gas with a median velocity of about 20 km s -1 . 
Interferometric CO images, H 2 O and the 18 km s -1 SiO 
masers, and dense-gas tracers such as thermal SiO emis¬ 
sion reveal a smaller (8" long) and younger (^ 200 year old) 
outflow along a northeast-southwest axis emerging from ra- 
dio source I orthogonal to the arc-minute-scale C O outflow 


(Beuther & Nissen 2008 Plambeck et al. 2009). The mo- 
mentum and kinetic energy content o f these flows is at le ast 

~ 1984) 


160 M e km s 1 and 4 x 10 4b ergs ( Snell et al.,, 

4 x 10 47 ergs ( Kwan fe Scoville||l976 ). Zapata etTaT 


to 


12009) 

presented a CO J = 2-1 interferometric study and found a 
dynamic age of about 500 years for the larger OMC1 out¬ 
flow. They noted that its impulsive nature and its structure 
is different from accretion-disk powered jets and collimated 
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protostellar outflows, and that it originated several arcsec- 
onds north of the OMC1 hot-core. 

Radio-frequency astrometry has shown that the three 
radio-emitting stars in OMC1, sources BN, I, and source 
n, have proper motions of 25, 13, and 26 km s -1 away 



its bipolar radio nebula which may shift source n’s emission 
centroid. 


Bally & Zinnecker (2005) proposed that the BN/KL ex- 
plosion may have been triggered by a massive star merger 
or the dynamical rearrangement of a nonhierarchical sys¬ 
tem of massive stars into a hierarchical system that resulted 
in their ejection from the OMC1 core as high-velocity, run¬ 
away stars. In this model, the disruption and ejection of cir- 
cumstellar disks and envelopes produced the BN/KL out¬ 
flow. The momentum and kinetic energy of the outflow and 
ejected stars came from the release of gravitational binding 
energy of a compact binary f ormed by the dynamic inter¬ 
action of three of more stars flPoveda et al. 1967[ Gualan- 


dr is et al. 2004 |Perets fe Subr||2012 Reipurth et al.pOlO 

Reipurth fe Mikkola||2012|). 

|Zapata et al. ( |2009p , Bally et al. ( 2011| ), and Goddi et al. 


( 201ip found supporting evidence for the dynamical rear- 
rangement model. In this scenario, the final N-body en¬ 
counter would have resulted in the formation of a compact, 
au-scale binary, most likely source I , and the ejection of 
radio source I and BN, and possibly source n. However, it 
is still possible that additional highly-embedded IR sources 
will be found in dense cores such as the 1 millimeter source, 


SMA1 (Beuther & Nissen 
tain massive stars. After a. 


2008). Some of these might con- 
17 radio source I has not yet been 


directly detected at any infrared wavelength. 

Over 30% of massive stars are ejecte d from their birth¬ 
places at high velocities (Gies & Bolton 1986). A large frac¬ 
tion of such runaway stars are thought to have been pro- 


duced by the dynamical ejection scenario (Poveda et al. 
1967). For example, a dynamical interaction i"s" thought to 


have ejected the classic runaway stars p Columbae and AE 
Aurigae with velocities of 108 km s -1 and 113 km s -1 , re- 
spectively, in opposite directions about 2.5 ±0.05 My r ago 


( Hoogerwerf et al. 12000 , 2001; Gualandris et al.|2004 ). This 
event left behind the colliding-wind X-ray binary, i Orio- 
nis, which is the brightest and most massive member of the 
5 Myr old NGC 1980 cluster located about 30' south of 


the Orion Nebula in the Orio n A molecular cloud (|Alves & 
Bouy||20l2 Bouy et ah]|2014 ). 


Proper motion measurements show that the fastest com¬ 
ponents in the OMC1 fingers have a dyna mic age of about 
500 years (Bally et al. 2011). However, Wu et al. (2014) 
used ALMA observations to argue for a somewhat older 


dynam ic age of <1,300 years for the BN/KL outflow. Bally 
et al. (2011) found that the H 2 and [Fe 11 ] finger proper 


motions trace-back to within a few arcseconds of J2000 
05:35:14.5, —05:22:23 locat ed between the current locations 
of radio sources I and BN. Gomez et al. (2008) found that 
in the Orion reference frame, the proper motions of radio 
sources I, BN, and n intersected within one arcsecond (< 


Article number, page 2 of 25 


400 au) of J2000 = 05:35:14.360, -05:22:28.70 about 500 
years ago, a few arcseconds from the apparent point of out¬ 
flow origin determined by tracing back the near-infrared 
proper motions of the fast ejecta. 


Tan (2004) proposed an alternative scenario in which 
the BN object was ejected about 4,000 years ago by a dy¬ 
namic interaction in the Trapezium cluster located in the 
center of the Orion Nebula. In this scenario, the OMC1 ex¬ 
plosion was triggered by the serendipitous passage of the 
BN object through the OMC1 core. Although requiring a 
highly unlike ly close encounter of BN with source I in the 


OMC1 core, Chatterjee & Tan (2012) show that the pa¬ 
rameters of the Trapezium and BN are compatible with 
this scenario. 

Though rare, the explosive outflow morphology of the 
OMC1 BN/KL outflow is not unique; other possible exam 
pies include G34.25+0.16 in the inner Galaxy flCyg; 

2008|), source G in W49 (|Smith et al.|f2009 

D 


et al. 


05506+2414 ( |Sahai et al.|2 008), 


anowski 

_ „ :009|), IRAS 

R21 (Zapata et al. 2013) 

urce of the molecular hy- 


and possibly the low-luminosity source 
drogen outflow MHO 2714 (GGD 34) in NGC 7129 (Eisloffel 


2000). However, Orion BN/KL is the nearest and least ob- 


scured, and thus most accessible for high-resolution studies. 
Here, we present 0.08 to 0.1 arc second resolution images 
of the entire OMC1 outflow complex in the 1.64 pm [Fe 11 ] 
and 2.12 pm H 2 narrow-band filters and a broad-band K s 
filter obtained with adaptive optics (AO) on the Gemini 
South 8 meter telescope. This data is combined with older 
AO-assisted observations obtained on Gemini North and 
natural seeing-limited images acquired with a variety of 
other telescope to re-measure proper motions in parts of 
the BN/KL outflow. 


2. Observations 

2.1. Gemini South GeMS 


The Gemini Multi-conjugate adaptive optics System 
(GeMS) at the Gemini South telescope on Cerro Pachon 
is the first and only sodium-ba sed multi-laser guide star 


(LGS) adaptive optics system flRigaut et al.| 2014 


2012 


Neichel et al. 2014 2013 Id’Orgeville et al. 2012|). GeMS 
works with a LGS constellation of 5-spots: 4 of the LGS 
spots are at the corners of a 60" square, with the 5-th po¬ 
sitioned in the center. The Adaptive Optics (AO) bench 
called Canopus is mounted on one of the f/16 Cassegrain 
ports. Gemini South Adaptive Optics Imager (GSAOI) is 
a wide-held 4096 by 4096 pixel (85" by 85" held of view) 
camera designed to work at the diffraction limit of the 8- 
meter telescope in the near-infrared. Three 85" diameter 
helds were observed in OMC1 between 30 December 2012 
and 28 February 2013 using GSAOI. Observations of each 
held were obtained through 1% bandpass narrow-band fil¬ 
ters centered on the 1.644 pm [Fell] and 2.122 pm H 2 emis¬ 
sion lines and the broad-band K s hlter. The corrected im¬ 
ages have FWHM diameters of 0.08" to 0.1", providing the 
highest angular resolution images of the BN/KL outhow 
ever obtained in the near-IR. 


Each held was imaged in each hlter with a 5 point dither 
pattern to hll-in gaps between the four 2048 by 2048 pixel 
arrays in GSAOI. Exposure times were 43.4 seconds per 
exposure for H 2 , 43.0 for [Fe II], and 15.0 for Ks. 10 ex¬ 
posures were taken in each hlter for a total of 430 seconds 
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on-source in the narrow-band filters and 150 seconds in the 
continuum filter. 

Data were processed with the Gemini GSAOI pipeline. 
However, additional astrometric corrections were requi red. 
Individual exposures were first registered to the |Muench| 


The two brightest [Fe II] bow shocks correspond to the 
Herbig-Haro objects HH 201 an d HH 210 located 60" north¬ 
west and 113" north of OMC1 ([Gull et al.||1973l |Mhnch fe 


et al. (2002) catalog sources to acquire a world coordinate 


51484 in 1999 (Kaifu et al 
Apac he Point Observato ry 3.5 meter on MJD 


2004 (Bally et al. 


53331 in 
The interval be- 


20111) were also used 
tween the Subaru and Gemini south observations is 4,839 
days. 


3. Results 

The 2013 epoch GSAOI images presented here reach the 
near-IR diffraction limit of an 8-meter telescope and provide 
the sharpest views obtained thus far of the entire OMC1 
BN/KL outflow. The combined H 2 , [Fe 11 ], and K s color 
image (Figure [l]) shows H 2 fingers tipped with [Fe 11 ] emis¬ 
sion extending from about 30" to 140" from the OMC1 
core. For the analysis of dynamic ages for various fea¬ 
tures, we assume that all features originated from J2000 
= 05:35:14.360, —05:22:28.70 (marked with a cyan cross in 
Figure [ 2 ]), the suspected location from which the BN o bject 
and radio source I were ejected about 500 years ago (Gomez| 
et ah| |2008 ) and within a few arc seconds of the suspected 
point of origin of the fingers determined from the intersec - 
tion point of the proper motion vectors (Bally et al.||2011). 
Figure [2] shows the H 2 image with the location of BN, I, 
n shown along with their proper motions. 


system with RMS pointing error ~ 0.1". A new catalog 
of relative star positions was generated from a preliminary 
aligned and co-added stack of images and used to derive a 
distortion map for GSAOI. The individual distortion cor¬ 
rected images were re-aligned and co-added to form the 
final mosaic in each filter. 


2.2. Gemini North Altair 

Between 2007 and 2009, the Gemini North 8 meter tele¬ 
scope was used to observe the OMC1 region on three oc¬ 
casions using the Altair AO system with the NIRI near-IR 
camera through 1% narrow-band [Fe II], H 2 , and broad¬ 
band K s filters. NIRI was used in a configuration which 
delivers a 40" field of view. In the narrow-band filters, a 
dithered set of five to ten 30 second exposures were ob¬ 
tained. A similar set of 10 second exposures were acquired 
in the K s filter. In 2007, only the ‘H 2 fingers’ region was 
observed as part of the commissioning of the Altair AO 
system using the NIRI camera. During 2008 and 2009, we 
intended to image a 3 x 4 point grid to cover the full ex¬ 
tent of the BN/KL outflow. However, only 5 and 8 fields 
were actually observed in 2008 and 2009, respectively. Only 
the ‘H 2 fingers’ field was observed during each of the three 
years. A summary of the observations is given in Table 1. 
The angular resolution of the NIRI images ranges from 0.1 
to 0.2". 

In the analysis presented here, the NIRI images were 
registered to the final GSAOI mosaic using IRAF tasks GE¬ 
OMAP and GEOTRAN applied to unsaturated field stars. 
Proper motions were determined by marking the photo¬ 
centers of features on the multi-epoch images. Images of 
the OMC1 outflow obtained with the_Subaru 8 meter tele¬ 
scope on MJD = 51484 in 1999 (jKaifu et al.|[2000), the 


Taylo r [1974;[Can to et al.|~L980| [A xon &Tay] [or 1984]). These 

shocks are visible on ground-based and Hubble Space Tele¬ 
scope images in [O i], Ha, [N Illand [S nj. Ho wever , they 
only exhibit faint H 2 emission ( (Graham et al.|[2003 ), indi¬ 
cating that they lie in the mostly atomic, photon-dominated 
region (PDR) located between the Orion Nebula’s ioniza¬ 
tion front and the background Orion A molecular cloud and 
OMC1 cloud core. The most prominent H 2 fingers consist 
of multiple [Fe 11 ] finger-tips trailed by H 2 wakes with an 
orientation of PA ~ 340° to 350°. This chain of shocks and 
wakes extends from about 50" to 135" from the ejection 
center. It consists of at least a dozen nested H 2 bow shocks 
tipped with [Fe 11 ] emission regions. The brightest [Fe 11 ] 
features are visible in the visual-wavelength emission lines 
commonly seen in Herbig-Haro objects such as [O 1 ], [S 11 ], 
[N 11 ], and Ha and have been designated as HH 205 through 
209. These HH objects are associated with the tips of the 
chain of H 2 wakes propagating toward position angle (PA) 
- 340° to 350°. 

More than 120 distinct wakes are visible in the 2.12 
/im H 2 images that exhibit near ly parallel walls and large 
proper motions along their axes ( Bally et al.|2011 ). Figure [ 3 ] 
shows a median filtered version of the 2013 GSAOI H 2 mo¬ 
saic created by convolving the image with a 51 pixel kernel 
(1"; each pixel is 0.02" on a side) using the IRAF function 
MEDIAN, and subtracting the result from the original im¬ 
age. Vectors were drawn from the suspected ejection site of 
radio sources BN and I (the coordinates are given above) to 
each H 2 or [Fe 11 ] fingertip. The dashed vector near the top 
marks a chain of H 2 knots and bow shocks with proper mo¬ 
tions nearly orthogonal to the northern fingers; this feature 
traces another flow originating east of the imaged field. This 
flow is also seen faintly at visual wavelengths. The H 2 emis¬ 
sion becomes too confused within ^30" of the suspected 
ejection location due to multiple overlapping features. The 
natura l seein g-limited 1999 epoch Subaru telescope image 
from Kaifu et al. (2000) and the 2005 epoch image from 
Bally et al. ( |201 1 ) were used to trace a few additional H 2 


an 


fingers beyond the boundaries of the GSAOI image in the 
southwest portion of Figure [3] 

In the northwestern part of the flow, the H 2 wakes range 
in diameter from 2" to 8" (7 x 10 14 to 3 x 10 15 cm) with limb- 
brightened rims less than 1" (< 3 x 10 14 cm) wide (Figure 
[4|. The half dozen major finger clusters in the northwest 
are up to 60" (^0.1 pc) long. The wakes in the inner part 
of the flow are narrower, tend to be shorter, and are more 
numerous resulting is a high degree of overlap along the line- 
of-sight. The H 2 emission tends to be fainter or disappears 
near the wake-tips where it is replaced by [Fe 11 ] emission. 

Most of the [Fe 11 ] emission in the BN/KL outflow orig¬ 
inates from the fingertips which are faint or invisible in H 2 . 
Dozens of fingertips are visible in the 1.64 /am [Fe 11 ] line. 
The two blue vectors in Figure [3] point to the two brightest 
[Fe 11 ] features which are also bright in visual wavelength 
shock-excited emission lines such as [O 1 ] and [S 11 ] and are 
designated as HH 201 and 210. These fingertips are very 
faint in H 2 . Figure [4] shows a color version of the ‘H 2 fin¬ 
gers’ region in 2013 from GSAOI showing H 2 (red), [Fe 11 ] 
(green), and the broadband light in K s (blue). 
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3.1. Proper Motions of Selected Knots 


Previous analyses of multi-epoch ground-based images have 
shown that the highest proper motion knots in the BN/KL 


outflow are most distant from OMC1 (Jones & Walker 1985 
Lee & Burton 2000; Bally et al. 2011). Figure [5] shows a dif¬ 


ference image ob tained by subtracti ng the 1999 epoch Sub¬ 
aru H 2 image of Kaifu et al. (2000) from our 2013 GSAOI 
H 2 image. Although the registration of the two images is 
not perfect due to small-scale distortions in the Subaru im¬ 
age (this image was assembled from a mosaic of smaller 
images), the overall pattern of expansion away from the 
OMC1 core is obvious. 


Comparison of the 2013 epoch images in the sub-field 
containing the brightest part of the PA ~ 340° to 350° 
finger with 2007, 2008, 2009 Gemini North NIRI images 
show large proper motions of the [Fe ll]-dominated finger¬ 
tips pointing away from the OMC1 core. This field contains 
several compact, high velocity H 2 knots with sub-arcsecond 
to arcsecond diameters and proper motions similar to the 
[Fe il] fingertips. Additionally, the H 2 wakes are expanding 
at right angles to a line connecting the fingertips to OMC1. 
The fastest transverse proper motions are found close to 
the fingertips and decline along the parts of the wake closer 
to OMC1. 


The most reliable proper motions were measured in this 
field using the AO-assisted images from 2013 and 2007 that 
are separated by 2,125 days. Figure [4] shows a color compos¬ 
ite image of this field. Figures [6] and [7] show difference im¬ 
ages in H 2 and [Fe il] formed by subtracting a de-distorted, 
intensity matched and registered 2013 GSAOI images from 
the corresponding 2007 NIRI images. Figures [8] and [9] show 
closeup views of the central part of the field shown in Fig¬ 
ures [6] and [7] that contains a bright bow shock. These images 
show 200 to 300 km s -1 proper motions of the fingertips in 
both H 2 and [Fe il]. 

Figure [8] shows that in addition to the large proper mo¬ 
tions of the fingertips, the H 2 wakes are spreading with 
velocities ranging from less than 20 to about 80 km s -1 in 
a direction orthogonal to the motions of the fingertips. The 
two regions marked W1 and W2 show the locations of slow 
and fast transverse spreading, respectively. The amount of 
transverse motion between the two epochs at these two lo¬ 
cations are indicated by the closely spaced dashed lines. The 
spreading velocity decreases with increasing distance from 
the fingertips. Figure [9] demonstrates that these spreading 
H 2 regions are located behind a fast-moving [Fe il] bullet 
whose locations in 2007 and 2013 are indicated. 

The closeup H 2 images of this region (Figures [6] and 
[8]) show a high-velocity compact clump (HVCC) which ex- 
nibits a proper motion of Vhvcc = 300 ± 20 km s -1 over 
the 6 year (1.85 x 10 8 second) interval between the acqui¬ 
sition of the Gemini North and South images in 2007 and 
2013. This HVCC is located at projected distance of 100" 
from the current location of source I and d e j ec u on = 97" 
(0.195 pc) from the suspected location of sources I and BN 
more than 500 years ago prior to their dynamic ejection. 
Assuming no deceleration, the dynamic age of this knot is 
tdyn = d eject ion, /V hvcc « 640 ± 30 years (see Table 2). 
Comparison of HST images in [O i] and [S il] taken on 
MJD 50170 show that the HVCC is located about 5" south 
of HH 207 which is also bright in our [Fe il] images. In 
[Fe il], HH 207 has a proper motion of 229 km s -1 (Ta¬ 
ble 2). Measured on Hubble Space Telescope images in the 
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visual wavelength range, HH 207 had proper motions of 
167 km s -1 in [S ill, 203 km s -1 in Ha, and 291 km s -1 in 
[N n] (D oi et al.| 2002). Using the velocity difference between 
the H 2 and [Fe nj proper motions of the HVCC and the 
[Fe il] emission associated with HH 207 (^68 km s -1 ) and 
the current projected separation of 5" indicates that the 
HVCC may catch up to the moving [Fe il] emission region 
in about 140 years (assuming that they are not displaced 
along the line of sight). 

Figure [8] shows that this HVCC has undergone consider¬ 
able transverse spreading over the past 6 years. The HVCC 
is either subject to photometric variations or experiencing 
significant “sideways splashing” of post-shock material. The 
morphological evolution of the HVCC suggests that it may 
have been decelerated significantly from its ejection veloc¬ 
ity. Several less-prominent HVCCs exhibit motions of 200 
- 300 km s -1 . Measured motions of a sample (subset) of 
fingertips, H 2 wakes, and compact, fast knots are given in 
Table 2. 

While the fingertips with 100 to 300 km s -1 motion 
(best seen in [Fe il]) are moving along a direction which 
points directly away from the OMC1 core, the H 2 wakes 
behind each [Fe il] feature are spreading at right angles to 
the fingertip proper motions. The wake spreading-velocities 
tend to decrease with increasing distance behind each fin¬ 
gertip. Within a few arcseconds of the fingertips, some H 2 
wakes show transverse expansion as fast as 80 km s -1 or¬ 
thogonal to the wake orientation and proper motion of the 
fingertip. Spreading speeds decline to under 20 km s -1 (the 
measurement limit) at distances of more than 10" to 15" 
behind the fingertip. 

The transverse spread is consistent with a simple model 
of a disturbance created by the passage of a high den¬ 
sity compact object through the medium as shown below 
by numerical simulations. Let y be the distance from the 
tip of a finger to the location where the wake half-width 
(width of the wake divided by 2) and transverse expansion 
speed are measured. The ratio of the transverse spreading- 
velocity, V P erp(y ) (one-half of the velocity with which the 
two sides of a wake are moving apart) to the fingertip 
proper motion, Upm, is comparable to the half-width of 
the wake at a given location, X per p{y ), divided by the dis¬ 
tance of that location from the fingertip, YpM{y )• Thus, 
Vperp{y)/Vp M « X perp (y)/Y PM (y). Both the fingertip and 
spreading velocities are highly supersonic. 

As shown by previous ground-based measurements, the 
largest proper motions are seen at the largest projected 
distance from the OMC1 core and the suspected locations 
of the BN object and radio source I about 500 years ago. 
The northwestern fingertips located at the greatest distance 
from the OMC1 core exhibit the largest proper motions 
(200 to over 300 km s -1 ) and tend to show a pattern of in¬ 
creasing maximum velocity with increasing distance from 
the core. However the region also contains a number of 
slower-moving features with proper motions ranging from 
under 100 to about 200 km s -1 . It is likely that these fea¬ 
tures have experienced significant deceleration as they in¬ 
teract with the surrounding medium. The fastest knots and 
fingertips have dynamic ages consistent with ejection be¬ 
tween 450 and 600 years ago. 


Doi et al.| ([2002) measured visual wavelength proper 
motions for Herbig-Haro (HH) objects protruding into the 
photon-dominated region (PDR) behind the Orion Nebula. 
Two of the brightest [Fe il] features correspond to HH 201 
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northwest of OMC1 and 210 north of OMC1. IDoi et ahl 
(2002) measured proper motions of 312 to 315 km s -1 for 


various components of HH 201 and and 309 to 425 km s 


for knots in HH 210. G rosso et al.| (|2006) detected X-rays 
from the wake of HH 210, which is the highest proper mo¬ 
tion finger in the OMC1 outflow and one of the relatively 
few visible at visual wavelengths, thereby demonstrating 
that at least some of the fingers contain hot, X-ray emit¬ 
ting plasma. Additional HH objects are associated with the 
Orion fingers including HH 205 to 209, and HH 601 to 
607 with proper motjpn velocities ranging from 100 to 300 
km s -1 (Doi et al. 2002). All of these features have large 
negative radial velocities . For example, H H 201 has Vlsr 
~ —260 to —284 km s -1 (Doi et al. 2004). 


3.2. Constraints on the ejection mechanism 

A number of models have been proposed to explain the 
[Fe il] fingers and H 2 wakes in the BN/KL outflow. The 
fingertips contain [Fe 11 ] knots and bow shocks with diam¬ 
eters of 1" 400 au) or less, estimated by measuring the 

separation between points where tangent lines on the lead¬ 
ing edges of the shock on either side of the tip are at an 
angle of 90° with respect to each other. This dimension is 
used as an estimator for the size-scale of the material pow- 
erin g each shock which m ay be the working surface of a 
jet ( Smith fe Rosen 2007|), a w ind which has exp erienced 
instabilities Q Stone et al.||1995), or a dense bullet (Allen &| 
Burton 1993). The tangent lines become increasingly par- 
allel (greater than 90°) in the wake. Most of these [Fe 11 ] 
bullets and the H 2 HVCCs are located more than 100" from 
their ejection site. 

The individual fingers resemble the shocks produced 
by accretion-powered jets or collimated flows powered by 
young stellar objects as th ey sweep up the ambient medium 
(Reipurth & Bally [2001). However the presence of many 
dozens of individual fingers and over 100 shocks with simi¬ 
lar dynamic ages moving away from the OMC1 core in all 
directions suggests that a single, brief event is more likely. 

The multitude of fingers and wakes could be powered by 
a fast wide-angle wind which experienced instabilities and 

t"" 


broke in to a m ultit ud e of p rotrusions (Stone et al 
McCaughrean & Mac Low|l997). Rayleigh-Taylor 


1995 


w 


stabilities can produce fingers of fast ejecta surrounded by 
slower clumps of dense gas if the wind velocity increases 
with time on a time-scale shorter than the crossing time 
of the wind in the wind-dominated region. Alternatively, a 
steady wind running into into a stationary medium with a 
density profile decreasing faster than r ~ 2 can be subject to 
RT instabilities. Neither an accelerating wind, nor a wind 
running into a rapidly decreasing density gradient is ex¬ 
pected to produce the approximately Hubble-flow type be¬ 
havior with the maximum velocity of the ejecta increasing 
linearly with increasing projected distance from the source. 

The fingers and wakes may be driven by high-density, 
com pact ejecta ( ‘bull ets’ or HVCCs) as originally suggested 


by Allen & Burton (1993). Similar bullets have been de¬ 
tected in other astrophysical contexts such as supernova 
remnants where hydrodynamic and cooling instabilities 


fragment an expanding ste llar envelope ( Fesen et al.||2011 
Milisavljevic fc Fesen||2Q13 ). 

Assuming the ejection occurred ^ 500 years ago at ap- 
proximately the location of the BN/KL infrared nebula 
(Bally & Zinnecker 2005 Zapata et al. 2009 Bally et al.| 


2011 Goddi et al. 2011), the most distant knots to the 
northwest provide the most stringent constraints on the 
ejecta properties. These knots have traveled ^ 0.28 pc in 
projection in ^500 years. Thus, the time-averaged velocities 
must be greater than ^550 km s -1 . For an explosive origin, 
the fingertip proper motions should increase linearly with 
increasing distance from the launch region. While a few H 2 
features have proper motions between 300 and 350 km s _1 , 
most show lower velocities, especially closer to the OMC1 
core. Because the faster motions are at least 20% slower 
than expected for a 500 year-old explosion, the ejecta have 
either decelerated or were ejected earlier. The absence of 
faster ejecta in [Fe 11 ] and H 2 could be a consequence of 
deceleration, an older time of ejection, or a selection effect 
resulting from excitation conditions in the shocks. [Fe 11 ] 
can be in excited into higher ionization stages and H 2 may 
be dissociated above a critical shock speed. Such ultra-fast 
material should still be revealed by their H 2 wakes. 

If the bullets originated between radio sources I and BN, 
they moved between 10 2 to 10 3 times their current diame¬ 
ters to their present location. The HVCC shown in Figure 
[8] had a diameter about about 0.2" (100 au) in 2007 and 
traversed a distance of about 100", about 500 times its own 
diameter, from the suspected ejection site between sources 
I and BN. Assuming that they are moving ballistically and 
not actively powered by jets or winds, momentum conser¬ 
vation requires that such bullets must be ^500 times denser 
than the environment. Such bullets are unlikely to contain 
sufficient mass to be bound by their self-gravity. They may 
be confined by ram pressure, otherwise they are likely to be 
expanding with a speed comparable to their internal sound 
speeds. 

If the bullets were formed from the disruption of circum- 
stellar disks envelopes within tens of au of the massive stars 
ejected from the OMC1 cloud core, their initial densities 
are expected to be n(H 2 ) ~ 10 10 to 10 15 cm -3 , compara¬ 
ble to the densities of disks and inner envelopes. Following 
dynamic ejection, bullets are likely to expand with a veloc¬ 
ity comparable to their an internal sound speed c s . In the 
absence of ram-pressure confinement by the medium the 
bullet is likely to sweep out a conical region given by the 
Mach angle, M ~ 2 c s /Vb where Vb is the bullet velocity. 
Using the observed sizes of the compact [Fe 11 ] and H 2 knots 
(40 to 400 au) and a distance greater than 0.2 pc from the 
ejection site implies c s < 0.6 km s -1 consistent with an H 2 
temperature of ~ 10 2 K . Although expansion will tend to 
cool the bullets below the surrounding medium, the intense 
IR radiation from OMC1 and heating by shock radiation 
will tend to keep them warm. It is likely that many bul¬ 
lets are moving into regions already shocked, heated, and 
partially evacuated by faster, previously ejected material. 
Such clumps will interact with a post-shock medium having 
a lower-density than the ambient medium, making it more 
likely for the bullets to be in free-expansion rather than 
being confined by ram pressure. 

The properties of the H 2 wakes constrain the density 
of the medium into which the suspected HVCCs are mov¬ 
ing. The wakes have widths that are an order-of-magnitude 
wider than the H 2 HVCCs and [Fe 11 ] fingertips (2" to 10" 
with most being near the lower-end of this range). Their 
formation requires that the heated layer behind the for¬ 
ward shock moving into the ambient medium has a cool¬ 
ing length la rger t han the clump (or the width of the jet 
beam) ( Blondin et al.|[l990 ). Figure [l0| shows a cartoon of 
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a dense bullet moving through a medium. A forward shock- 
heated ambient medium streams around the moving clump 
and drives lower-velocity shocks sideways into the ambient 
medium. A much slower, reverse shock propagates into the 
clump and compresses its leading edge. 

The post-shock temperature immediately behind a 
shock is given b y T vs =_3pV 2 /16k where V s is the shock 
speed. Blondin et al. (fl990) give the cooling distance as 
VTMr 


dcool = = 9fiV^jMn 0 A(T ps ) 


4.5 x 10 lu t?oo/«o 
where p is the mean molecular weight of the pre-shocked gas 
particles, A (T ps ) is the cooling function, Vioo is the shock 
velocity in units of 100 km s -1 , and no is the pre-shock 
particle density. More recent numerical studies of the cool¬ 
ing function give d coo i & 5.5 x 10 17 VAn 4 /no for 80 km s -1 
< V s < 1,200 km s _1 (Draine 2011). Thus, a V s = 300 
km s _1 shock moving into a density no = 10 4 cm -3 has 
a cooling length L coo i ~ 7 x 10 15 cm (470 au). Thus, for 
ambient medium densities between 10 3 cm -3 and 10 4 cm -3 
typical for the ‘Integral Shaped Filament’ in the Orion A 
molecular cloud located behind the Orion Nebula, these 
cooling lengths correspond to 1.1" to 11" at the distance of 
Orion, and are larger than the HVCCs and [Fe il] knots. 

For a dense HVCC driving a forward shock into the 
medium with speeds of 150 to 500 km s -1 , the hot (~1 
to 10 MK) plasma will splash sideways to produce a wide 
bow-shaped wake. In the OMC1 rest frame the 30 to 300 
year cooling time corresponds to a wake-length L wa j^ e ~ 
3 x 10 16 to 3 x 10 17 cm, the latter being comparable to the 
lengths of fingers in the north and west parts of the OMC1 
outflow. The sideway expansion of the hot plasma into the 
surrounding medium can then drive a slower ( V s id e <100 
km s -1 ) side shock. The H 2 emission is likely to be produced 
in the swept-up, compressed, and shock heated layer behind 
the shock. However, intense UV radiation produced by the 
shock hot plasma may also excite H 2 molecules outside the 
wake, resulting in fluorescent excitation. 

The reverse shock propagating back into a dense bullet 
will be much weaker than the forward shock owing to the 
bullet’s higher density. For a forward shock moving into a 
medium with density no with velocity Vs, the speed of the 
reverse shock moving into the bullet having a density np (in 
the frame of the bullet) is Vp ~ Vsirio/np) 1 ^ 2 . For a density 
contrast ratio np/n 0 ~ 10 3 , the reverse shock propagating 
back into the bullet will have speed Vp ~ 0.031(5 = 10 
km s -1 for a forward shock speed of 300 km s -1 . 


3.3. Numerical modeling 

We use the Eulerian adaptive mesh refinement (AMR) hy¬ 
drodynamics + N-body code, Enzo (Bryan et al. 2014) to 
model the propagation and evolution of bullets ana the 
shocks they produce as they migrate through t he_a m bient. 
The simulat ion set up is analogous to that of |Silvia et al. 


2012 ) in which an over-dense, spherical cloud is em- 
d in a uniform ambient medium. For the simulation 


(2010 

beddo 

presented here, we have assumed an adiabatic equation of 
state that ignores any effects radiative cooling may have 
in this scenario. More detailed calculations using a realistic 
cooling curve for the post-shock gas will be presented in a 
future paper. 

In order to track the largest possible head-tail struc¬ 
ture as the cloud is disrupted, without losing information 
due to flows that might exit the computational domain, we 
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give the ambient medium a uniform velocity such that it 
washes over the cloud, rather than moving the cloud itself. 
As the medium flows past the cloud, it drives shocks into the 
cloud interior and ablates material from the cloud-medium 
interface - producing a tail of material that is carried down¬ 
stream. In order to capture the evolution of the cloud over a 
large physical volume while minimizing the computational 
cost of the simulation, we use a coarse-resolution root grid 
with a physical cell size of ~78 au per cell and eight levels 
of adaptive mesh refinement such that the cloud structure 
is resolved at a peak resolution of ^0.3 au per cell. The full 
computational volume is 7500 x 2500 x 2500 au. 

Figure [TT shows a comparison between one of the H 2 fin¬ 
gertips and the numerical model showing the density struc¬ 
ture of a dense, supersonic cloud that has moved approx¬ 
imately 500 times its initial diameter through a uniform 
medium. At the start of the simulation, the medium is ini¬ 
tialized with a number density of n — 10 3 cm -3 , a temper¬ 
ature of T = 20 K, and uniform velocity of v = 250 km s -1 
directed along the x-axis of the simulation domain. The 
cloud is initialized with a density 10 4 times higher than 
that of the ambient medium, a temperature of 100 K, and 
a radius of r = 10 au. 


Figures [l2| through [I6| made using the yt toolkit (Turk 


et al. 201l)7sh ow the density, x, y,z -velocity structure, and 


temperature of the bullet after 100 years of evolution. For 
this simulation, the initial velocity of the cloud with the 
respect to the medium is too slow to reach the observed 
positions given an ejection 500 years ago, but is compatible 
with an earlier ejection time. Additionally, the simulation 
shows significant cloud deceleration and disruption in just 
120 years following ejection. Together, this indicates that 
the bullet over-densities are likely to be higher than 10 4 , 
or that the ISM densities are lower than used in the mod¬ 
els. Prior to the dynamic ejection event, the forming stars 
probably produced accretion-powered outflows which would 
have carved cavities in the surrounding cloud through with 
ejecta would be moving. Additionally, some bullets may be 
propagating through cavities and wakes created by the lead¬ 
ing edge of the ejected debris field. 

The inclusion of cooling in the post-shock regions will 
tend to lower the transverse expansion of the shock-heated 
plasma to produce a narrower wake and a lower sideways 
splashing velocity. Cooling of the shock-heated layer pro¬ 
duced by the slower, reverse shock moving back into the 
bullet will also tend to result in extensive fragmentation 
as portions of the cloud rapidly cool and condense - a phe^ 


nomenon observed in the cloud-crushing simulations of (Sil¬ 


via et al. 2010 2012). Such fragments could have consid¬ 


erably longer survival times than the structures present m 
the simulation presented here. Other physical phenomena 
not included in this numerical work (e.g. magnetic fields, 
thermal conduction, plasma viscosity) may also aid in cloud 
survival. The addition of such physics is left to future work. 


3.4. Other YSOs in the field 

The line-of-sight to the OMC1 BN/KL outflow contains 
both lower-mass protostars embedded in or near the OMC1 
cloud core, and more mature young stellar objects in the 
foreground Orion Nebula. Several dozen YSOs fall into the 
field-of-view presented here. Several of these YSOs drive 
outflows which can be seen in our images. 
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The most prominent outflow crosses the northern part 
of the BN/KL outflow around Declination —05:31 at posi¬ 
tion angle, PA ~ 240° to 250° and indicated by a dashed 
arrow near the top of Figure [3j Com parison w ith the [O i] 


and [S n] HST images presented by Doi et al. 


that this flow consists of a bow shock located at J2000 


(2002) show 


5:35:06.8, —5:18:27. Do i et al.j (2002) label this feature as 
115-101 and 115-102 in their Figure 4 where it is shown 
as part of HH 603 and measured a proper motion toward 
PA ~ 230° with a velocity of 70 to 75 km s -1 . In the H 2 
image, a larger but fainter bow shock is located north of 
115-101 and 115-102 at J2000 = 5:35:11.6, -5:20:54. In our 
H 2 images, this 2" radius bow shock shows a proper motion 
parallel to 115-101 and 115-102 with a velocity of about 70 
km s -1 . This feature is also faintly seen in the HST [O 1 ] 
and [S 11 ] images obtained in 1994. Assuming that the [O 1 ], 
[S 11 ], and H 2 emission arise from the same region, these 
images confirm the mostly west-southwest motion. Finally, 
the [O 1 ] image shows a very faint bow-shaped feature at 
J2000 = 5:35:08.9, —5:18:10 with a morphology and ori¬ 
entation consistent with the above proper motion vectors. 
Apparently, this set of three west-facing and west-moving 
shocks traces an outflow which originates from a yet-to-be 
identified YSO located northeast of the OMC1 core in the 
Integral Shaped Filament of dense gas in the Orion A cloud. 

The variable star V2270 Ori (05:35:15.394, 
—05:21:14.11) located 71" north of BN along the line-of- 
sight to the eastern edge of our images drives a bipolar 
[Fe 11 ] outflow at PA ~ 34° (Figure 17). The flow consists 
of a pair of bow-shaped features l.fr southwest and 2.1" 
northeast of the star. The southwest portion of the flow 
forms a jet extending several arc seconds beyond the 
southwest bow. A curved, less-collimated jet-like feature 
extends to the northeast. This morphology is reminiscent 
of a jet produced by coverging flows at the end o f a nozzle 
(Canto & Rodriguez 1980 Canto et al. |1981 ). Spectra 
were obtained at the Apache Point Observatory 3.5 meter 
telescope using the Triplespec spectrometer (these results 
will be presented fully in a forthcoming paper). These 
spectra show that the southwestern jet is redshifted with 
a mean radial velocity of about +30 km s -1 with respect 
to the rest frame of the Orion Nebula [Fe 11 ] emission. The 
northeast lobe of the V2270 [Fe 11 ] jet is blueshifted with a 
radial velocity of about +30 km s -1 . The velocity changes 
at the position of the star, indicating that V2270 Ori is 
the driver of this outflow. 

There is a silhouette disk system at 5:35:14.692 
—5:22:20.36 about 8" east of BN (Figure [l8|) . It is most 
evident in the [Fe 11 ] images, but also seen in the H 2 and 
Ks filters. The arcsecond-scale bipolar nebula perpendicular 
to the shadow likely traces reflected continuum light from 
the central star embedded in and shadowed by the nearly 
edge-on disk. The shadowing occurs at position angle 9° 
implying a disk axis at about 99°. 


4. Conclusions 

New images of the OMC1 BN/KL outflow, obtained with 
the Gemini South MCAO laser guide star adaptive optics 
system are presented. The images, obtained in the 2.12 fi m 
H 2 line, the 1.64 /im [Fe 11 ] line, and in a broad-band K s 
filter nearly reach the diffraction limit of the 8 meter tele¬ 
scope, about 0.08". Complete coverage of the outflow’s H 2 


wakes and [Fe 11 ] fingertips reveal over 100 individual shocks 
in unprecedented detail. 

These images are combined with 2007 to 2009 epoch 
images obtained with the Gemini North AO system, Altair, 
and the near-IR camera NIRI, and earlier epoch natural 
seeing limited images to measure new proper motions of 
selected features and to search for morphological changes 
in the BN/KL outflow. 

The observed morphologies of the fingers and wakes are 
compared with numerical hydrodynamic simulations of a 
compact, high density bullet moving through a lower den¬ 
sity ambient medium with a Mach number of about 10 3 . 
The main results of this investigation are: 

— Several dozen [Fe 11 ] bright fingertips and a few sub- 
arcsecond diameter H 2 knots show the highest proper 
motions. The proper motions are 150 to 300 km s -1 at 
the projected north end of the outflow complex. 

— The survival of H 2 in a knot (the HVCC) moving with 
a proper motion or order 100 km s -1 indicates that it 
must have a density much higher than the surrounding 
medium. Over the last 6 years, this feature has shown 
dramatic morphological changes as it approaches a high 
velocity [Fe 11 ] knot located several arcseconds father 
downstream. The H 2 emission from this high-velocity 
compact clump (HVCC) has spread laterally (orthogo¬ 
nal to its proper motion), expanding from a diameter 
less than 0.1" to more than 0.4", possibly indicating 
the rapid evolution of the reverse shock moving into the 
clump 

— Many of the wakes show measurable expansion with a 
velocity of up to ~ 80 km s -1 orthogonal to the proper 
motions of the fingertips which produce them. The com¬ 
pact structure of the [Fe 11 ] and H 2 emission of the fin¬ 
gertips suggests that they are powered by sub arcsecond 
to arcsecond scale bullets which have a much greater 
density than the medium into which they are moving. 
It is argued that the density in the bullets must be more 
than three orders-of-magnitude greater than the back¬ 
ground density. 

— A simple model of bullet propagation is presented 
in which a fast forward shock produced mega-Kelvin 
plasma whose cooling length is much longer than the 
diameter of the bullet. The expansion of this plasma 
drives shocks into the medium at right angles to the 
bullet motion to produce < 100 km s -1 shock which 
can sweep up the ambient medium without completely 
destroying H 2 . These sideways shocks driven by the ex¬ 
pansion of the shock-heated medium are responsible for 
the intense H 2 emission in the OMC1 BN/KL outflow. 
The reverse shock moving into the bullets is slow and 
may excite some detectable H 2 emission. 

— Preliminary numerical simulations with the Enzo code 
using an adiabatic equation of state reproduce the ob¬ 
served morphology and kinematics of the fingers. The 
numerical simulations reproduce the observed sideways 
splash, the spreading of the post-shock material orthog¬ 
onal to the propagation of the bullets, and the fragmen¬ 
tation of the leading edge of the high-velocity debris. 
Further numerical modeling is needed to fine-tune the 
model parameters to better match the densities and dy¬ 
namical ages of the observed OMC1 shocks. 
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Fig. 1 A wide-field image showing the OMC1 outflow in H 2 (orange) and [Fe ii](cyan). 
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Fig. 2 A wide-field image showing the OMC1 outflow in H 2 with proper motions of the BN object, radio source I, and 
radio source n superimposed. The lengths of the red-solid vectors are proportional to the motions measured by Gomez 
et al. (2008) with the the lengths of the vectors in arcseconds equal to the motion in km s -1 (e.g. 10" corresponds to a 
motion of 10 km s _1 ). The ejection center as determined by radio proper motions is shown by a cross. The locations of 
the disk shadow shown in Figure [l8| and the source of the [Fe 11 ] jet shown in Figure pT| V22T0 Ori, are indicated. 


Article number, page 10 of 


25 





Bally et al: Orion GeMS GSAOI 


Dec. (J2000) 

30.0 23:00.0 30.0 22:00.0 30.0 21:00.0 30.0 -5:20:00.0 


J3 

> 



o 

o 



Fig. 3 The outer system of fingers (vectors) in the 2.122 /i m H 2 GSAOI image. The dashed vector shows the background 
flow whose motions are towards the west. The blue vectors mark the direction from the explosion center to the two [Fe 11 ]- 
dominated shocks, HH 201 and HH 210. The background image has been median filtered as described in the text. 
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? ig. 4 A color image showing the K s image (blue), the 1.644 p m [Fell] image (green), and 2.122 pm H 2 (red) in the ‘H 2 
ingers’ field using the GSAOI 2013 data. 
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Fig. 5 A 2.122 /im H 2 difference image showing proper motions in the field imaged by Kaifu et a. (1999) with the Subaru 
telescope in 1999. The images shows the difference between 2013 and 1999 epoch data. Residual distortion corrections 
result in imperfect registration. 
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Fig. 6 A 2.122 /am H 2 image showing proper motions in the ‘H 2 fingers’ field. The image shows the difference between 
images obtained in 2013 with Gemini S using GSAOI and in 2007 with Gemini N using NIRI. 
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Fig. 7 A 1.644 /im [Fell] image showing proper motions in the ‘H 2 fingers’ field. The image shows the difference between 
images obtained in 2013 with Gemini S using GSAOI and in 2007 with Gemini N using NIRI. 
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Fig. 8 A closeup view of the main H 2 finger and the HVCC in the ‘H 2 fingers’ field showing the difference between images 
obtained in 2013 with Gemini S using GSAOI and in 2007 with Gemini N using NIRI. 
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Fig. 9 Same as Figure [7] but for [Fe il] . 
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Fig. 10 A cartoon showing the shock structure formed by a hypersonic, dense bullet moving through a cold medium. 
The flow vectors are shown in the rest-frame of the bullet. The forward shock velocity Vo is nearly identical to the bullet 
velocity and heats the post-shock ambient medium to at temperature ^ 10 6 K. The reverse shock propagated back into 
the bullet with a speed Vr ~ Vb^o/n#) 1 / 2 , where no is the density of the ambient medium and tir is the bullet density. 
The transverse expansion of the hot, high-pressure, post-shock medium sweeps-up the colder, lower-pressure ambient 
medium to form the H 2 wake. Only the transverse component of the vector V S id e is shown. The fluid velocity of the 
material swept-up by the expansion of the hot, shocked ambient medium is the vector sum of the pre-shock ambient 
medium (as seen in the reference frame of the bullet), and the velocity of the sideways splash produced by the expansing 
hot plasma. 
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Fig. 11 A numerical simulation of a spherical object with a density 10 3 times higher than the background density, moving 
with a velocity 10 3 times faster than the sound speed in the medium. The frame shows the projected density distribution 
after the objects has moved about 10 3 times its initial diameter. The simulation projection is superimposed next to one 
of the fingers shown in Figure 1. 
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Fig. 12 A density-weighted projection of gas density from the Enzo simulation described in Section 3.3 at t = 100 years. 
Projected values are computed such that p pro j = J2 p! / Pi where pi is the density, i denotes the cell number, and p 

i i 

is the density in cell i. 
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Fig. 13 A slice through the shock model showing the x-velocity from the Enzo simulation described in Section 3.3 at t = 
100 years. The velocities correspond to motions along the line-of-sight. Owing to the nature of the simulation, since the 
bullet is initialized at rest the medium is flowing past it, small x-velocities signify material that has not been significantly 
accelerated. Or, when frame-shifted to the case of a moving bullet, material that has not been significantly de-accelerated. 
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Fig. 14 A slice through the shock model showing the ^/-velocity from the Enzo simulation described in Section 3.3 at t 
= 100 years. The velocities correspond to motions along the line-of-sight. 
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Fig. 15 A slice through the shock model showing the z-velocity from the Enzo simulation described in Section 3.3 at t 
= 100 years. The velocities correspond to motions perpendicular to the motion of the bullet and the line-of-sight. 
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Fig. 16 A density-weighted projection of gas temperature from the Enzo simulation described in Section 3.3 at t = 100 
years. Projected values are computed such that T pro j = ^2[Ti x pi] j ^ pi, where Ti is the projected quantity, pi is the 

weighting quantity, and i denotes the cell number. As discussed in the text, actual post-shock temperature in slices that 
are not weighted by density reach temperatures in excess of 10 6 K. 



Fig. 17 The V2270 [Fe il] jet (blue) and H 2 emission in the region (orange). 
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Fig. 18 The near-IR disk shadow located 8" east of BN. [Fe il] emission is shown in blue; H 2 emission is shown in orange. 
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Table 1 Near-IR Observations of Orion OMC1 


Year 

a(J2000) 

(5(J2000) 

Field Size 

MJD 

Comments 

2007 

05:35:12.8 

-05:20:53 

50" 

x 50" at 45 0 

54165 

H 2 fingers; Gemini N 

2008 

05:35:11.9 

-05:20:47 

50" 

x 50" at 45 0 

54753 

H 2 fingers; Gemini N 

2009 

05:35:11.9 

-05:20:47 

50" 

x 50" at 45 0 

55138 

H 2 fingers; Gemini N 

2013 

05:35:13.2 

-05:20:37 

84" 

x 84" at 0° 

56290 

OMCl N; Gemini S 

2013 

05:35:12.5 

-05:22:04 

92" 

x 91" at 0° 

56323 

OMC1 BN; Gemini S 

2013 

05:35:17.5 

-05:22:48 

92" 

x 91" at 0° 

56323 

OMCl SE; Gemini S 


Note. — The H 2 , [Fe 11 ], and K s images were taken on the same night. 
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Table 2 Selected OMC1 bullets and bow shocks 


Feature 

a(J2000) 

(05^ 35 m +) 

<5(J2000) 
(-5° + ) 

AX 

n 

V 

(km s -1 ) 

PA 

(°) 

D 

(") 

^ ~dyn 

(yrs) 

Comments 

H 2 1 

12.625 

20:55.22 

0.88 

297 

346 

97 

643 

HVCC H 2 knot 

2 

12.058 

20:37.68 

0.43 

145 

324 

116 

1577 

Ahead of H 2 1 

3 

12.270 

20:37.82 

0.84 

283 

344 

114 

794 

11 

4 

12.351 

20:38.84 

0.54 

182 

350 

112 

1213 

11 

5 

12.965 

21:16.37 

0.33 

111 

338 

74 

1314 

Cluster NE of H 2 1 

6 

12.989 

21:17.73 

0.53 

179 

347 

73 

801 

11 

SI 

11.532 

21:02.24 

0.30 

101 

233 

- 

- 

Background flow 

S2 

11.616 

20:53.94 

0.29 

98 

227 

- 

- 

ii 

S3 

11.614 

20:54.45 

0.54 

182 

218 

- 

- 

ii 

W1 

12.543 

20:52.03 

0.240 

80 

- 

- 

29 

Spreading H 2 wake 

W3 

12.642 

20:56.45 

0.138 

46 

- 

- 

100 

11 

[Fe II] 1 

12.521 

20:50.12 

0.68 

229 

345 

102 

863 

HH 207; Fe II] bow 

2 

12.165 

20:39.09 

0.81 

273 

345 

116 

835 


3 

12.274 

20:37.96 

0.79 

267 

354 

114 

839 


4 

12.353 

20:38.94 

0.55 

186 

349 

112 

1184 


5 

12.898 

20:34.12 

0.59 

199 

1 

116 

1146 


6 

12.932 

20:34.87 

0.61 

206 

3 

115 

1097 


7 

13.066 

20:35.39 

0.61 

206 

355 

114 

1088 


8 

13.052 

20:39.01 

0.61 

206 

354 

112 

1069 


9 

13.010 

20:40.81 

0.56 

189 

350 

97 

1009 


10 

13.308 

20:54.78 

0.41 

138 

354 

96 

1367 


11 

13.386 

20:56.16 

0.53 

179 

353 

95 

1043 


12 

13.354 

20:53.14 

0.33 

111 

357 

98 

1736 


13 

12.592 

20:56.23 

0.27 

91 

354 

97 

2096 


14 

12.795 

21:05.50 

0.49 

165 

354 

87 

1036 


15 

11.502 

21:07.20 

0.92 

310 

338 

94 

596 


16 

11.680 

21:10.86 

0.71 

240 

326 

89 

729 


17 

11.792 

21:12.81 

0.41 

138 

329 

86 

1225 



Note. — [1]: The projected distance from the suspected point of ejection of the stars BN and source I about 500 years ago, D, 
is measured from J2000 = 05:35:14.35, —5:22:28.5 (Gomez et al. 2008). 
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